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Abstract 

Zn and Ni impurities in the hole-doped high-temperature superconductors 
are known to have strong effects on thermodynamic and transport properties. 
A recent scanning tunneling microscope study of Zn-doped Bi2212 (Pan et al, 
Nature 403, 746 (2000)) has provided high-resolution images of the local den- 
sity of states around non-magnetic impurities in dx2_y2-wave superconductors. 
These pictures contain detailed information about the spinor wave functions 
u{r) and v{r) of bound states with energy Eq ~ A/30, centered at the Zn- 
sites. We show that this type of wave function follows from the solutions of 
the Bogoliubov-de Gennes equations for dx2_y2-wave superconductors. 
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Introduction: Impurity-doping of anisotropic superconductors has turned out to be a 
valuable probe, setting them apart from conventional superconductors with an s-wave order 
parameter. For example, doping of high-temperature superconductors, such as YBCO, 
LSCO, and BCCO, with non-magnetic impurities has contributed greatly to establish the 
underlying dx2_y2-wave order parameter in these systems. [|T]-|^ In particular, an analysis of 
the thermodynamic and transport properties in these compounds suggests that Zn-impurities 
can be modeled with a scattering potential in the unitary limit, On the other hand, little 
is known experimentally about the local real-space structure of the impurity bound states 
around the Zn-sites, in spite of numerous theoretical studies on this question. 0-^] Very 
recently. Pan et al. have provided detailed scanning tunneling microscope (STM) images 
of the local tunneling conductivity, fixed at ±£'0? where ~ A/30 is the binding energy 
of the impurity bound state around a Zn site in Bi2212. The corresponding wave function 
exhibits a fourfold symmetry, associated with the underlying CUO2 lattice. 

In the context of conventional s-wave superconductors, this type of bound state, in this 
case associated with a magnetic impurity, has been predicted by Shiba. |Q Recently, this 
idea has found experimental support by STM studies in the vicinity of Fe and Gd impurity 
atoms on the surface of Nb superconductors. 



The objective of this paper is to present a simple solution of the Bogoliubov-de Gennes 
equations for dx2_y2-wave superconductors |lTT| , p^ in the presence of a single impurity. This 



corresponds to an experimental situation where the impurity concentration is dilute enough, 
such that interactions between the impurities can be neglected. Within this approach, 
we obtain a wave function which resembles very accurately the fourfold symmetry which 
has recently been observed around Zn-impurities in Bi2212 p[, where the energies of the 
associated bound states are very small {\Eq\//S. <^ 1). We believe that this image of the 
local tunneling conductivity around the Zn-site, as seen by STM, provides strong evidence 
for dx2_y2-wave symmetry. Therefore, this method can also be used to test anisotropic 
order parameters in organic superconductors, and may serve as a straightforward tool to 
explore the underlying symmetry in other new superconductors where the nature of the 



order parameter is still under debate. 

We also study the case of Ni impurities which may be modeled by a weaker scattering 
potential, close to the Born limit. |13| In this case, we find that although |f (r)p is not much 
different from the Zn-case, |u(r)p is rotated by 7r/4 with respect to the strong-scattering 
limit corresponding to Zn. Finally, we find that for impurities with an impurity potential of 
intermediate strength, the lowest harmonic in |'u(r)p is suppressed, and therefore a dominant 
eight-fold symmetry may be expected instead of the four-fold symmetry, observed in the weak 
and strong scattering limits. 

Bogoliubov-de Gennes equations and bound state wave functions: The Bogoliubov-de 
Gennes (BdG) equations for dx2_y2-wave superconductors in the continuum limit are given 

by 00 

Eu{r) = (^-^ - - y(r)) u{r) + ^^A{dl - dl)v{r), (1) 

Ev{r) = - (^-^ - /i - Vir)^ v{r) + ^A{dl - dl)u{r\ (2) 

where /i is the chemical potential, and V{r) > is the impurity potential, centered at the 
site r = 0. 

Let us first consider the case of Zn impurities where the energy of the bound states is 
known to be very small [§], Eq ~ 0. In the following, we use a variational ansatz for the 
solutions of the BdG equations: 

u{v) = A exp (-7r) [Mppr) + V2(3Ji{pFr) cos (40)) , (3) 
v{y) = V2Aaexp {—^r)J2{pFr) cos (20), (4) 

where Ji{z) are Bessel functions of the first kind, and pp is the Fermi momentum. [A is the 
global normalization factor for the wave functions which can be neglected.) Note that for 
Bi2212 it is believed that pp ~ 0.7A~1. This implies that ppa ~ 2. [|l4l The coefficients a, 
/3, and 7 are determined variationally. Inserting Eqs. 3 and 4 into Eqs. 1 and 2, we find 

E = K -V - ^Aa, 



Ea = -Ka-^Ail + ^), (5) 
E(3 = K(3- ^Aa, 

where 



K 



drr [dj-exp {—'yr)Ji{ppr))'^ + (/exp {—^r)Ji{pFr)/r)^ 



2m /q°° (irr (exp (— 7r) J/(pi?r))^ ' mpp 



fx ^ (6) 



In general, K in Eq. 6 depends on which Bessel function Ji{pf'^) is used. However, in the 
hmit 7/pf <^ 1, -ft' reduces to '-/^/{mpF) for all Ji{pFr) with / <^ Pf/i- In the approximation 
for V, only the dominant s-wave component (/ = 0) of the scattering potential V^(r) has 
been considered. In the usual convention, which is used in this paper, a positive sign of 
V{r) corresponds to the attraction of electrons by the impurity. In impurity-doped Bi2212 
the charge carriers are holes, and hence the opposite sign has to be chosen. 

For the case of a Zn- impurity (strong scattering limit), we may assume that ~ 0. p| 
This gives K ^ V/2, V ~ 3A/2 ^ V^A, a ~ V2{V - K)/A, and (3 ~ -^2. These choices 
yield the approximate parameter set a ~ 1, /3 ~ — 1/-\/2, and 7 ~ PpiV^PpO^^^^ — O.SpF- 
The tunneling conductance is given by ||15[ 



dl , , ,2 feV - Eo\ 2 , 1 2 feV + EQ^ ^ 



^(r) (xsech^ l^^j I^Wr + sech^ l^^j l^^"^' " 

Thus, at small temperatures, the local tunneling conductance around the impurity site is 
dominated by |'u(r)p for a fixed binding energy Eq and by |f (r)p for —Eq. In Figs. 1 and 2 
images of |'u(r)p and |f (r)p are shown. Both |M(r)p and |f (r)p have a four-fold symmetry, 
and extend in the directions of the Cu-0 bonds. This strongly resembles the images seen by 
the STM experiments. ^ Furthermore, weaker higher-harmonic satellite peaks are observed 
in |u(r)p, along the (±7r/4, ±7r/4) directions, in accordance with the experiments 0. Note 
that with our sign convention for the scattering potential ^(r), the roles of |M(r)p and |f (r)p 
are interchanged for the hole-doped cuprates. Thus, in Bi2212 our |M(r)p corresponds to a 
bound state at —Eq and |f (r)p to a bound state at Eq. 



Now let us turn to the case of Ni impurities. Ni is considered to be a weak scatterer. 
Therefore we may assume that V{r) —>■ 0, and consequently a ~ 2(3{K — E)/A and (3 ~ 
72(1 - 4(^2 _ j^2^ /A2)-l. So, if > _^ AV4 the coefficient (3 turns out to be positive, 
which is most likely the case for Ni. Indeed, a sufficiently large magnitude of a implies that 
^2 > ^ and E ~ V3A/2. This gives a ~ {3 ~ 1/^2, and 7 ~ O.lpF- The 

squares of the wave functions ■u(r) and f (r) are shown in Figs. 3 and 4. Comparing with 
the corresponding images for the strong scattering limit, we note that |'y(r)p is qualitatively 
similar in both cases, whereas |M(r)p appears to be rotated by 7r/4 with respect to the case 
of Zn impurities. 

For impurities with a scattering potential ^(r) of intermediate strength, the coefficient 
j3 in the wave function u{y) is found to be greatly enhanced. Consequently, the higher- 
harmonic term in u{v) (Eq. 3), containing a cos (40)-modulation dominates the tunneling 
response at positive binding energies +£'0. In this case, an eight-fold instead of a four- fold 
symmetry should be observed in |M(r)p, as shown in Fig. 5. On the other hand, the shape of 
|f (r)p is not affected. The weaker satellite peaks in Fig. 1 can be viewed as strong-coupling 
precursors of this phenomenon. It would be interesting to test this particular prediction 
of an eight-fold symmetry in |M(r)p experimentally, using appropriate candidate impurity 
atoms with intermediate scattering strengths. 

Conclusions: In summary, we observe that (i) using the Bogoliubov-de Gennes equations 
for single impurities in dx2_y2-wave superconductors in the continuum limit is a very useful 
approach, (ii) Both Zn and Ni impurities produce bound states with a four-fold symmetry, 
localized around the impurity sites, (iii) Impurities with intermediate scattering strength 
may have a bound state with a dominant eight-fold symmetry. The simple analysis of the 
BdG equations we have presented in this work provides a semi- quantitative picture of the 
spatial structure of bound states in anisotropic superconductors. 

Recent experiments suggest that the electron-doped high-temperature superconductors 
NCCO and PCCO also have a d-wave order parameter. [|T^-|T^ On the other hand, the 
effects of Zn and Ni doping in the electron-doped systems appear to be opposite from the 
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hole-doped case, i.e. Ni impurities lead to a stronger suppression of the superconducting 
ordering temperature than Zn doping. |19[ It is therefore likely that the roles of Zn and Ni 
are interchanged in the electron-doped high-Tc cuprates. Furthermore, there are indications 



that the k — (ET)2-salts have a d-wave superconducting order parameter as well. p0|j21 



Therefore, similar studies of impurity bound states in these compounds would be of great 
interest . 

We thank Hyekyung Won for useful discussions. S. H. acknowledges the Zumberge 
foundation for financial support. 
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-4 -2 2 4 
FIG. 1. Spatial variation of the local tunneling conductance, centered at a strong-scattering 
impurity, such as Zn, in a dx2_y2-wave superconductor. In this figure, the dominant contribution 
[^(r)^ at the positive bound state resonant frequency -{-Eq is shown. 



-4-2 2 4 
FIG. 2. Spatial variation of the local tunneling conductance, localized around a 
strong-scattering impurity, such as Zn, in a dx2_y2-wave superconductor. Here the dominant con- 
tribution |v(r)p at the negative bound state resonant frequency —Eq is shown. 
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-4 -2 2 4 
FIG. 3. Spatial variation of the local tunneling conductance, centered at a weak-scattering 
impurity, such as Ni, in a dx2_y2-wave superconductor. In this figure, the dominant contribution 
|'u(r)p at the positive bound state resonant frequency +Eq is shown. 
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FIG. 4. Spatial variation of the local tunneling conductance, localized around a 
weak-scattering impurity, such as Ni, in a dx2_y2-wave superconductor. Here the dominant contri- 
bution |f(r)p at the negative bound state resonant frequency —Eq is shown. 
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FIG. 5. Spatial variation of the local tunneling conductance, centered at an impurity of 
intermediate scattering strength. Here the dominant contribution [^(r)^ at the positive bound 
state resonant frequency +Eq is shown. In contrast to the weak and strong scattering limits, a 
dominant eight-fold symmetry is observed. 
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